Bond-Stretching-Phonon Anomalies in Stripe-Ordered Lai 69Sro.3iNi04 
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We report a neutron scattering study of bond-stretching phonons in Lai.69Sro.3iNi04, a doped 
antiferromagnet in which the added holes order in diagonal stripes at 45° to the Ni-0 bonds. For the 
highest-energy longitudinal optical mode along the bonds, a softening of 20% is observed between the 
Brillouin zone center and zone boundary. At 45° to the bonds, a splitting of the same magnitude is 
found across much of the zone. Surprisingly, the charge-ordering wave vector plays no apparent role 
in the anomalous dispersions. The implications for related anomalies in the cuprates are discussed. 
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There is resurgent interest in the role of phonons with 
respect to the high-temperature superconductivity found 
in layered copper-oxides m. Particularly striking are 
the anomalies in high-energy optical modes observed 
by neutron scattering in La2_a;Sra;Cu04 |2| ||, Q| and 
YBa2Cu306-i-£!; 1^, H ^, 1^. There have been various 
speculations as to vifhether the observed phonon anoma- 
lies might be related to instantaneous charge inhomo- 
geneities, particularly those in the form of stripes M. Re- 
cently, detailed analyses of electron-phonon interactions 
in a dimerized stripe phase have been reported |9[ p^ . 

One way to learn about the effect of charge stripes 
on lattice dynamics is to study a model system with 
well-defined stripe order. Here we present the first 
single-crystal study, to our knowledge, of the bond- 
stretching phonon modes in such a system, specifically 
La2-2;Sr2;Ni04 with x — 0.31. The stripe order in Sr- 
doped nickelates has been characterized in detail by neu- 
tron diffraction, and the most recent summary of re- 
sults is given in |ll|. For x > 0.22, the crystal struc- 
ture is tetragonal, consisting of a body-centered stack- 
ing of Ni02 planes. Within the Ni02 planes, the charge 
stripes run diagonally along either [110] or [110] di- 
rections, at 45° to the Ni-0 nearest-neighbor bonds, 
which extend along [100] and [010] directions. (In con- 
trast, the charge stripes observed in superconducting 
Lai.g_a;Ndo.4Sr2;Cu04 run parallel to the Cu-0 bonds 
[p2|; however, for x < 0.06, the stripes inferred to exist 
in La2-xSra;Cu04 have the diagonal orientation of the 
nickelates [ [l3| , |l4|.) The maximum transition tempera- 
tures for charge-stripe and magnetic order occur at a; = g 
|ll| , |5|. Well below the charge-ordering temperature, 
Tco, the nickelates have very large resistivities [1^, [l7[ , 
consistent with all of the added holes being localized in 
charge stripes. 

We focus on the highest-energy longitudinal optical 
modes propagating along the [100] and [110] directions. 
These modes show only weak dispersion in stoichiomet- 
ric La2Ni04 [18|. In contrast, we observe a softening 
of 20% along [100] on moving from the Brillouin zone 



center to the zone boundary, quite similar to that found 
in La2-2:Sr:rCu04 [g[ Q|; along [110], a splitting of the 
same magnitude is observed over much of the zone and, 
in particular, at the zone boundary. These results are 
important for two reasons. 1) The observed anomalies 
must be associated with the local charge inhomogeneity. 
They are induced by the hole doping, and the holes are 
localized in the stripes. 2) There is no evidence that the 
charge-ordering wave vector, q^.^, plays a special role. If 
the phonon anomalies were related to collective phase 
fluctuations of the charge stripes, as in a conventional 
charge-density- wave system [n9|, then one might expect 
them to appear at qco- With the absence of a collec- 
tive signature, it seems likely that the dominant effects 
involve local interactions between charge and lattice fluc- 
tuations. 

Our La2_a;Sra;Ni04 crystal, grown by the floating-zone 
method, is cylindrical, with a diameter of 6 mm and 
length of 30 mm; the Sr concentration of x = 0.31 was 
confirmed by inductively-coupled plasma analysis (with 
an uncertainty of ±0.01). At room temperature, the lat- 
tice parameters of the tetragonal unit cell are a = 3.84 A 
and c — 12.70 A. The charge and spin ordering tran- 
sitions were confirmed to be consistent with previous 
work[jril (magnetic and charge-order transitions of ap- 
proximately 160 K and 235 K, respectively) by neutron 
diffraction measurements performed at the JRR-3M re- 
actor in Tokai, Japan. The inelastic-neutron-scattering 
measurements of phonons were performed on the IT 
triple-axis spectrometer at the Orphee reactor of the Lab- 
oratoire Leon Brillouin in Saclay, France. The spectrom- 
eter is equipped with a Cu (111) monochromator and a 
pyrolytic graphite (PG) (002) analyzer, each of which is 
both vertically and horizontally focused. The analyzer 
was set to detect neutrons with a final frequency of 3.55 
THz (14.7 meV). A PG filter was placed after the sam- 
ple to minimize unwanted neutrons at harmonic wave- 
lengths. 

Figure 1(a) shows scans of scattered intensity versus 
excitation frequency for phonons propagating in the [100] 
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FIG. 1: Neutron scattering measurements of phonons in 
Lai.69Sro.3iNi04 at T = 10 K. (a) Scattered intensity vs. 
frequency for reduced wave vectors, q = (ft, 0,0), varying 
from Brillouin zone center {h = 0) to the zone boundary 
{h — 1). Curves through the data are fits to gaussian peak 
shapes on top of a constant background. Scans have been off- 
set vertically for clarity, (b) Similar results for q = {h,h,0); 
dashed lines indicate the constant background, (c) and (d) in- 
dicate the positions in reciprocal space where the scans (iden- 
tified by letters) of (a) and (b), respectively, were measured; 
dashed lines indicate Brillouin zone boundaries. Measure- 
ments were performed in various equivalent zones, with vary- 
ing L-component, in order to avoid spurious peaks due to 
accidental Bragg scattering by the sample. 



direction, parallel to the in-plane Ni-0 bonds. The scans 
are taken at wave vectors along the line q = {h, 0, 0) 
in the first Brillouin zone, moving from the zone center 
{h — 0, bottom) to the zone boundary {h — 1, top), 
vifhere the components of the wave vector are measured 
in reciprocal lattice units, (27r/a, 27r/a, 27r/c). (The effec- 
tive zone boundary for a single Ni02 plane is at /i = 0.5; 
the fact that the actual boundary is at /i = 1 results 
from the body-centered stacking of the layers.) Each 
scan is dominated by a single, well-defined peak, that 
disperses from approximately 21 THz at zone center to 
less than 18 THz half-way across the zone, and then back 
up again. The lines through the data points are fitted 
gaussian peaks on top of a background that is taken to 
be independent of wave vector and frequency. 

Figure 1(b) shows scans measured along the [110] di- 
rection. Because of twinning of the stripe domains, we 
simultaneously probe phonons propagating parallel and 



perpendicular to the stripes. Again, we focus on the 
mode that starts at about 21 THz at zone center; the 
mode at 14 THz involves bond-bending motion, and does 
not exhibit any notable doping-dependent behavior. [The 
intensity variation seen in the figure for this mode is as- 
sociated with its sensitivity to the L-component of the 
wave vector Q, which varies from point to point as noted 
in Fig. 1(d).] On moving from zone center to zone bound- 
ary, we note that the highest-frequency mode varies little 
in frequency, but it loses intensity. Concommitant with 
this, signal appears in the 17-18 THz regime, growing in 
strength as the zone boundary is approached. 

Using the gaussian fits as a smoothed version of the 
data, we present our results as intensity plots in Fig. 2(a) 
and (b). Here, the intensity has been multiplied by I'/Qf, 
(Q|| is the component of Q along the phonon propaga- 
tion direction) to correct for sensitivities of the neutron- 
scattering cross section. Dramatic differences are ob- 
served in the dispersions of the stripe-ordered system 
compared to those in stoichiometric La2Ni04, which are 
indicated by the gray circles llq]. Contrary to naive 
expectations based on the Peierls-distortion model [ p^ , 
there is no particular anomaly at the charge-ordering 
wave vector qco = (0.31,0.31,0) [see Fig. 2(b)]. In- 
stead, we find a splitting of the mode along [110] over 
a substantial part of the zone. The size of the disper- 
sion in the [100] direction is about the same as the [110] 
splitting, and is essentially identical to the magnitude 
of the dopant-induced softening reported for supercon- 
ducting Lai.85Sro.i5Cu04 [g|,|4|. The softened modes are 
consistent with a previous study of the phonon density- 
of-states in La2-2:Sr^Ni04 ^^, where a dopant-induced 
peak was observed at ~ 75 meV (~ 18 THz). 

Integrating the corrected intensity data of Fig. 2(a) 
and (b) over frequency yields the open circles shown in 
Fig. 2(c) and (d), respectively. The experimental re- 
sults are compared with intensities calculated from the 
interatomic-potential model used to describe the phonons 
in La1.gNiO3.93 [Q. (An overall normalization of the cal- 
culated to the measured intensities has been applied in 
each panel.) The general consistency between the mea- 
sured and calculated integrated intensities suggests that 
the distribution of weight for the anomalously-softened 
bond-stretching modes does not extend significantly be- 
low 16 THz. 

To interpret our results, let us consider some sim- 
ple phenomenological models. The high-energy bond- 
stretching phonons involve the motion of the in-plane 
oxygens between the much heavier nickel atoms. Sup- 
pose we assume the nickel atoms to be infinitely heavy 
and consider only nearest-neighbor forces; then each oxy- 
gen acts like an Einstein oscillator, with no dispersion 
of its vibrational frequency. Correcting for the finite 
mass of the nickel would result in a slight decrease in 
frequency towards the zone boundary; however, account- 
ing for Coulomb repulsion between neighboring oxygen 
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FIG. 2: (color) False-color images of intensity as a function 
of frequency and q for (a) q — {h, 0, 0) and (b) q = (ft, h, 0). 
The plotted intensity corresponds to the fitted curves (minus 
the background) from Fig. 2. qco indicates the approximate 
charge-ordering wave vector. The gray circles indicate the 
phonon dispersions measured in stoichiometric La2Ni04|J18[. 
(c) Intensity vs. q integrated over frequency along the [100] 
direction. Open circles are data, stars are intensities calcu- 
lated from the interatomic-potential model for La1.9NiO3.g3 
[pi. (d) Similar plot for the [110] direction. 



ions would counter that with an increase in frequency 
at the zone boundary, where neighboring oxygens move 
opposite to one another. Ignoring these corrections, a 
simple Einstein model gives a rough approximation of 
the measured dispersions in stoichiometric La2Ni04 [111] 
[see Fig. 2(a) and (b)]. 

To go further, we will restrict ourselves to a lin- 
ear Ni-0 chain model. The softening of phonons at 
the zone boundary can be described by introducing 
a force between nearest-neighbor O ions with a nega- 
tive force constant, /cbr HI, [see Fig. |3|(a)]. This 
force phenomenologically incorporates a particular type 
of electron-phonon coupling. An example of the disper- 
sion from such a model is shown in Fig. |3|(c). The eigen- 
vector for the softened zone-boundary mode, sketched in 
Fig. ||(e) , involves linear breathing motion of the O about 
the Ni. 

Next, suppose that fcbr acts only between one pair of 
oxygens out of three [see Fig. ||(b)]. The result is a split- 
ting of the modes, as illustrated in Fig. ^(d). Note that 
the intensity of the softened mode is strong at the zone 
boundary, but goes to zero at zone center. This effect is 
achieved without a dimerization of the lattice. As shown 
in Fig. 0(f), the eigenvector of the softened mode at zone 
boundary involves an isolated, local breathing motion of 
a pair of O about a Ni site. The behavior of this mode 
does not require coherence between motions in different 
unit cells. 

It seems likely that both of these models have some 
degree of relevance to our observations; however, it is 
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FIG. 3: (color) (a) Sketch of a force-constant model for a 
linear Ni-O array with a nearest-neigbor force constant of 
strength k, and an effective force constant fcbr between next- 
nearest-neighbor oxygens, (b) Model with fcbr acting only 
between one out of every three oxygen pairs, (c) Calculated 
intensity vs. q = {2n/a)h and frequency for the model in (a) 
with {2k/moy^^/2Tr = 21 THz and (-2fcbr/mo)^/V27r = 
9 THz. (d) Calculated intensity for model (b) with 
(2fc/mo)^/V27r = 21 THz and {-k^^JmoY^'^ /2-k = 9 THz. 
(e) Sketch of the eigenmode ai h — 0.5 for model (a), corre- 
sponding to a linear breathing mode, (f) Softened eigenmode 
at /i = 0.5 for model (b), corresponding to a local breathing 
mode. 



difficult to associate them in a consistent or unique way 
with specific details of the measured dispersions. Of 
greater interest is the nature of the electron-phonon 
coupling modelled by fcbr- In model calculations for 
La2._a;Sra;Cu04, Falter and Hoffmann |2|] have empha- 
sized the importance of ionic charge fluctuations in re- 
sponse to atomic displacements and the irrelevance of 
Fermi-surface nesting effects of the type once proposed 
by Weber |2^. Given the experimental evidence for 
poor electronic screening of phonons in optimally-doped 
cuprates pJ , extended Coulomb interactions should also 
be important. 

In the present case of the nickelate, little metallic 
screening is expected due to the gap of > 0.1 eV in the op- 
tical conductivity at low temperature [|l^, ^ . The bond- 
stretching modes involve polar fluctuations of the nega- 
tive oxygen ions against the positive nickel ions, and so 
should couple to charge fluctuations. The lowest-energy 
channel for charge fluctuations must involve the dopant- 
induced holes, and since the holes are segregated into 
stripes, charge fluctuations must be associated with some 
form of stripe fluctuations. 

In a conventional charge-density-wave (CDW) system, 
one might expect the mode along [110] to couple to fluc- 
tuations of the phase of the CDW with respect to the lat- 
tice liy]. In the present case, the absence of a character- 
istic wave vector suggests that local charge fluctuations 
may be more relevant for the bond-stretching modes. Lo- 
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FIG. 4: Comparison of scans for q = {h, h, 0) at temperatures 
of 10 K (filled circles) and 250 K (open circles). The lines are 
the same as in Fig. 2. The upper scans have been shifted 
vertically by 200 counts. 



cal fluctuations were predicted by Yi et al. ^^ in a 
random-phase-approximation treatment of Hartree-Fock 
stripes. Nevertheless, our results do not rule out the pos- 
sibility of weak, low-energy, collective phase fluctuations 
of the stripes. 

Do the effects we observe depend on static ordering 
of the stripes? Figure ^ shows a comparison of phonons 
at zone-center and zone boundary along [110] measured 
at 10 K and 250 K, which is above the charge-ordering 
temperature. One can see that at 250 K, the zone-center 
mode has softened slightly but remains strong and well- 
defined; however, at the zone boundary, the intensity re- 
mains spread out in frequency. Thus, the anomalies do 
not disappear with the static order. 

Rather similar phonon softenings were observed in the 
sample identified as La1.9NiO3.93 ||l^, |2^. Despite the 
specified stoichiometry, that sample exhibited superlat- 
tice peaks of the type (7J/2, K/2, L/2) with H, K, and L 
odd, consistent with the stage-2 ordering of interstitials 
found in La2Ni04.io5 P7| . The latter compound exhibits 
evidence of incipient stripe correlations, but no stripe or- 
der |2^. This is further evidence that local interactions 
are sufficient to yield the phonon anomalies. 

To summarize, we have observed doping-induced 
anomalies in bond-stretching modes of Lai.69Sro.3iNi04, 
a compound in which the doped holes order in stripes 
at low temperature. The anomalies persist over a sub- 
stantial portion of the Brillouin zone, with no obvious 
signature of collective stripe fluctuations. The splitting 
of the mode along [110] is similar to recent observations in 
superconducting Bao.eKg. 46103 [g^, where charge inho- 
mogeneity is also expected to be relevant. These results 
provide support for a connection between the phonon 



anomalies observed in the cuprates and instantaneous 
charge inhomogeneity. 
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